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We study the dynamical exponent z for the directed percolation depinning (DPD) class of models 
for surface roughening in the presence of quenched disorder. We argue that z for (d+ 1) dimensions 
is equal to the exponent dmin characterizing the shortest path between two sites in an isotropic 
percolation cluster in d dimensions. To test the argument, we perform simulations and calculate z 
for DPD, and dmin for percolation, from d = 1 to d = 6. 
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Recently the growth of rough interfaces has been the 
object of many theoretical, numerical, and experimental 
studies, fueled by the broad interdisciplinary aspects of 
the subject M. Applications are as diverse as imbibi- 
tion in porous media, fluid-fluid displacement, bacterial 
colony growth, fire front motion, and the motion of flux 
lines in superconductors R] . 

In general, a d-dimensional self-afline interface, de- 
scribed by a single-valued function h(x,t), evolves in a 
(d -|- l)-dimensional medium. Some form of disorder, 77, 
affects the motion of the interface leading to its roughen- 
ing. Two main classes of disorder have been discussed in 
the literature. The first, called "annealed", depends only 
on time. The second, "quenched" disorder, is frozen in 
the medium. 

Continuum equations, such as the Kardar-Parisi- 
Zhang (KPZ) equation [g[, have been remarkably suc- 
cessful in describing roughening for the case of annealed 
disorder ||l|. For the quenched disorder case, several mod- 
els were proposed with a view toward explaining exper- 
imental results for which the roughness exponent a is 
significantly larger than the predictions for annealed dis- 
order (for reviews, see e.g. ^). Here W ^ L°', where W 
is the interface width and L is the system size. 

For one class of models, the static properties of the 
interface in (1 + 1) dimensions are obtained exactly by a 
mapping, at the depinning transition, onto directed per- 
colation (DP) [g|j]. In higher dimensions the mapping is 
to directed surfaces (DS) B — for (1-f 1) dimensions, DP 
and DS are equivalent. This class of models is referred 
to as directed percolation depinning (DPD) [pj. 

A recent numerical study ||7|, confirmed by analytical 
results [pj , showed that the DPD class of models can be 
described by a stochastic differential equation of the KPZ 
type with quenched disorder 
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/i(x, t)=F + uV'^h + \{Vhf + ?7(x, h). 
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Here F is the driving force and 77 (x, h) is the quenched 
disorder acting to roughen the interface; above the depin- 
ning transition, F > Fc leads to an interface that moves 
indefinitely with a constant velocity. Reference |0] finds 



for the DPD class of models that the coefficient of the 
nonlinear term. A, diverges at the depinning transition 

Because of the presence of a diverging nonlinear coef- 
ficient in (nh, the application of the functional renormal- 
ization group to the calculation of the exponents has not 
been possible ||l|. For this reason, the only existing esti- 
mates of the exponents are from simulations, or from the 
mapping of the scaling properties of the pinned interface 
to DP (d = 1) or DS (d > 1) |-|. Since no mapping 
of the dynamics of Eq. (^ or of the DPD models has 
been found, no theoretical estimates for the exponents 
characterizing the dynamics of the roughening process 
has been made. One such exponent is the dynamical ex- 
ponent z, which characterizes the lateral propagation of 
the interface perturbations, as is defined by pO[ |, 

t-rl (2) 

where t is the time needed for a perturbation to spread 
over a longitudinal distance rn . 

In this Letter we argue that z for the DPD universality 
class can be identified in {d+l) dimensions with the per- 
colation exponent dmin for the shortest path for isotropic 
percolation in d dimensions nw. We support this rela- 
tion by numerical calculations of both z and dmin, up to 
dimension (6-1- 1). Our work implies an upper critical 
dimension for the dynamics, dc + 1 = 7, above which the 
mean field result, z = 2, becomes exact. 

In the DPD model the growth occurs on a discrete lat- 
tice and the disorder is modeled by considering that each 
cell has a probability p of being "blocked" ^-^ . Since 
the model was developed to study imbibition, we will re- 
fer to the growing, invading, region as "wet" , and to the 
remaining region as "dry" . At time i = 0, we wet all cells 
at the bottom of the lattice. Then, at each time step, we 
wet all dry unblocked cells that are nearest-neighbors to 
a wet cell. To retain a single-valued interface, we im- 
pose the auxiliary rule that all dry blocked cells below a 
wet cell become wet as well. These cells we call "eroded 
blocked cells" , and this procedure is referred to as erosion 
of overhangs [||-|| . If the concentration of blocked cells 



is small, the interface propagates forever. As the con- 
centration of blocked cells increases, large portions of the 
interface become pinned by fragments of DP clusters in 
d = 1 or by fragments of DS for d > 1. The characteristic 
longitudinal dimension of these fragments, fii , diverges as 
p approaches the critical threshold pc- When <^|| becomes 
comparable with the system size L, the interface eventu- 
ally becomes completely pinned by a spanning DP path 
or DS. Just below pc, almost all of the interface is pinned 
except for a few unblocked points which move along the 
interface creating new sites for growth. We address the 
behavior of the system only in its critical state, i.e. on 
length scales smaller than ^|| , where ^|| ~ (pc — p)""" • 

In order to find z, we study how perturbations caused 
by a single unblocked cell propagate over the interface. 
At each time step a certain set of cells become wet by 
incoming fluid, caused by a single unblocked cell at time 
t = 0. In analogy with invasion percolation, we call this 
set of cells the percolation shell. For each time step i, we 
compute the average radius of gyration of the percolation 
shell r{t). Since t is the time needed for a perturbation 
to spread over a distance ry, t obeys Eq. (g). 

For the case d = 1, all shells are confined between 
the old directed path that spans the system at i = 
and a new pinning path that will block the growth after 
some time. The region between these two paths is ef- 
fectively one-dimensional, since the vertical distance be- 
tween them scales as the perpendicular correlation length 
of DP, ^±_ ^ [pc — p)~^^ . Hence v± < i/^\ implies 
^_l/C|| ^ as p ^ Pc- For any cell on the interface that 
becomes wet at time t, one can find the cell from which it 
was invaded at the previous time step, and recreate the 
sequence of invasion events that leads from the initial cell 
to any given cell on the interface (Fig.nl). The trajectory 
of this sequence follows the upper pinning path and is 
effectively one dimensional. Its length i scales as its av- 
erage end-to-end distance ry . On the other hand, £ is 
equal to the time t needed to reach the end of the path. 



Hence t 



and we conclude from 



that z = 1. This 



conclusion is supported by our simulations (Table ffl) . 

For the case d > 1, we must consider the region 
bounded by two self-afhne, single- valued, DS (Fig.||) [|j. 
This region is effectively d-dimensional, since ^±/^\\ -^ 0. 
Hence, the shortest path leading from the initial point 
to any point of this region is effectively confined to a d- 
dimensional horizontal hyperplane. This shortest path 
has to avoid blocked cells in this hyperplane, as does 
the shortest path of isotropic percolation. For isotropic 
percolation it is known that the length of the shortest 
path £ scales with the Euclidean end-to-end distance r as 
£ ~ r'^raia^ -pjj^g similarity between the geometrical prop- 
erties of the paths in DPD and isotropic percolation leads 
us to propose 
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We present our results for the exponents z and d,„iii in 
Table |. 

It is well known that for isotropic percolation the up- 
per critical dimension is dc = 6, i.e. for d > dc the mean 
field result, dmin = 2, becomes exact [|ll|. This suggests 
an upper critical dimension, dc -I- 1 = 7, for the dynamics 
of the DPD models which are in the universality class of 
Eq. (0), and that z = 2 for d + 1 > 7. 

Since the dynamics of Eq. ffl) and the models in the 
DPD universality class are connected to isotropic perco- 
lation, while the static properties are mapped to DP or 
DS, it is possible that the upper critical dimension deter- 
mined in this study may be valid only for the dynamics. 
In fact, it is possible that dc for the static properties 
may not exist. Suppose, e.g., that a one-dimensional 
object, such as a self-avoiding walk, is embedded in a 
d-dimensional space. We expect that as d is increased 
the interactions between the different parts of the object 
decrease. At a certain d = dc, these interactions can be 
neglected, and the exponents become those of the ideal 
non-interacting case. In contrast, when the dimension 
of the object is not fixed but increases with d, as in the 
case of DS in which the object is one dimension smaller 
than the space, we expect to move away from the non- 
interacting limit. In fact, the analytical solution of the 
DPD model in the Cayley tree suggests that the upper 
critical dimension for the statics might be oo pi. 

In summary, we present an argument that identifies the 
dynamical exponent z for the DPD universality class with 
the fractal dimension of the shortest path in isotropic 
percolation, dmin. This result leads us to identify the di- 
mension (6-1-1) as the upper critical dimension for the 
dynamics of the DPD universality class. 
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To test the argument leading to (0), we performed sim- 
ulations for both DPD and percolation for d = 1 to d = 6. 
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FIG. 1. Illustration of the dynamics of the DPD model for 
(1 + 1) dimensions, (a) Schematic representation of a region 
defined by two pinning paths. The heavy circle indicates the 
origin for the invasion, the thin arcs represent the positions 
of the invading front at successive times, and the dashed line 
represents schematically the path for the invasion, (b) Simu- 
lation results for invasion after 2^" time steps starting from 
a single cell near the center. We show the invaded region at 
a sequence of times which are multiples of 128. Regions in- 
vaded at later times are displayed in darker shades of gray. 
The path from the origin to the latest invaded point is shown 
in black. Although this path displays some fluctuations in the 
vertical direction, they can be disregarded since i^n > i/±, so 
as p — > Pc C±/C|| -^ 0. Thus the distance propagated by the 
invading front is proportional to time. Since ix ~ ^, we can 
conclude that z — dmin ~ 1. 



FIG. 2. Illustration of the dynamics of the DPD model 
for (2 -I- 1) dimensions, (a) Schematic representation of the 
xy projection of the region defined by two pinning self-affine 
DS. The heavy circle indicates the origin for the invasion, the 
thin arcs represent the xy projections of the invading front at 
successive times, and the dashed line represents schematically 
the path for the invasion, (b) Simulation results for invasion 
after 2^" time steps starting from a single cell located to the 
left of the center. We show the xy projection of the invaded 
region at a sequence of times which are multiples of 128. Re- 
gions invaded at later times are displayed in darker shades 
of gray. It is visually apparent that it takes a long time to 
invade some regions close to the origin because the path to 
that position (shown in black) appears to be a fractal curve of 
dimension greater than one. The fluctuations in the vertical 
direction can be disregarded since we know that £,±/^\\ —> 0. 
We find that the path can be identified with the shortest path 
(the "chemical distance" ) of isotropic percolation, and that its 
length scales with the linear distance r to the point as r"^""". 



FIG. 3. (a) Scaling with time of the horizontal length of 
a DPD cluster in (d-f 1) dimensions grown from a single cell. 
Shown is a double logarithmic plot of time t as a function 
of r|| , which is the average of the parallel components of the 
radius of gyration of the shell. The asymptotic slope is z. 
(b) Double logarithmic plot of the shortest path £ in isotropic 
percolation as a function of the Euclidean distance r. The 
asymptotic slope is dmin. Note that, after some transient be- 
havior, a transition to a power law scaling occurs. For higher 
dimensions, the power law scaling is affected by finite-size 
effects for larger times. 



TABLE I. Dynamical exponent z for the DPD model in 
id + 1) dimensions and the shortest path exponent dmin for 
isotropic percolation for a d-dimensional cubic lattice of L'* 
sites. The results indicated by an asterisk are exact, while 
the remaining values were calculated in our simulations by 
the study of the consecutive slopes of the linear regime in 
Fig. 3. At the critical dimension dc = 6, one should not 
expect to find the exact result dmin = 2 because logarithmic 
corrections are generally present. The system sizes used in 
the simulations range from L = 4096, for d = 2, to L = 16, 
for d = 6. 
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